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Abstract 

Gravitational wave solutions to Einstein's equations and their generation are examined in D- 
dimensional flat spacetimes. First the plane wave solutions are analyzed; then the wave generation 
is studied with the solution for the metric tensor being obtained with the help of retarded D- 
dimensional Green's function. Due to the difficulties in handling the wave tails in odd dimensions 
we concentrate our study in even dimensions. We compute the metric quantities in the wave zone in 
terms of the energy momentum tensor at retarded time. Some special cases of interest are studied: 
first the slow motion approximation, where the D-dimensional quadrupole formula is deduced. 
Within the quadrupole approximation, we consider two cases of interest, a particle in circular orbit 
and a particle falling radially into a higher dimensional Schwarzschild black hole. Then we turn 
our attention to the gravitational radiation emitted during collisions lasting zero seconds, i.e., hard 
collisions. We compute the gravitational energy radiated during the collision of two point particles, 
in terms of a cutoff frequency. In the case in which at least one of the particles is a black hole, we 
argue this cutoff frequency should be close to the lowest gravitational quasinormal frequency. In this 
context, we compute the scalar quasinormal frequencies of higher dimensional Schwarzschild black 
holes. Finally, as an interesting new application of this formalism, we compute the gravitational 
energy release during the quantum process of black hole pair creation. These results might be 
important in light of the recent proposal that there may exist extra dimensions in the Universe, 
one consequence of which may be black hole creation at the Large Hadron Collider at CERN. 
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I. INTRODUCTION 



One expects to finally detect gravitational waves in the forthcoming years. If this hap- 
pens, and if the observed waveforms match the predicted templates, General Relativity will 
have pass a crucial test. Moreover, if one manages to cleanly separate gravitational wave- 
forms, we will open a new and exciting window to the Universe, a window from which one 
can look directly into the heart of matter, as gravitational waves arc weakly scattered by 
matter. A big effort has been spent in the last years trying to build gravitational wave 
detectors, and a new era will begin with gravitational wave astronomy [1, 2]. What makes 
gravitational wave astronomy attractive, the weakness with which gravitational waves are 
scattered by matter, is also the major source of technical difficulties when trying to develop 
an apparatus which interacts with them. Nevertheless, some of these highly non-trivial 
technical difficulties have been surmounted, and we have detectors already operating [3-5]. 
Another effort is being dedicated by theoreticians trying to obtain accurate templates for 
the various physical processes that may give rise to the waves impinging on the detector. We 
now have a well established theory of wave generation and propagation, which started with 
Einstein and his quadrupole formula. The quadrupole formula expresses the energy lost to 
gravitational waves by a system moving at low velocities, in terms of its energy content. The 
quadrupole formalism is the most famous example of slow motion techniques to compute 
wave generation. All these techniques break Einstein's equations non-linearity by imposing 
a power series in some small quantity and keeping only the lowest or the lowest few order 
terms. The quadrupole formalism starts from a flat background and expands the relevant 
quantities in i?/A, where R is the size of source and A the wavelength of waves. Perturbation 
formalisms on the other hand, start from some non-radiative background, whose metric is 
known exactly, for example the Schwarzschild metric, and expand in deviations from that 
background metric. For a catalog of the various methods and their description we refer the 
reader to the review works by Thorne [6] and Damour [7] . The necessity to develop all such 
methods was driven of course by the lack of exact radiative solutions to Einstein's equations 
(although there are some worthy exceptions, like the C- metric [8]), and by the fact that 
even nowadays solving the full set of Einstein's equations numerically is a monumental task, 
and has been done only for the more tractable physical situations. All the existing methods 
seem to agree with each other when it comes down to the computation of waveforms and 
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energies radiated during physical situations, and also agree with the few available results 
from a fully numerical evolution of Einstein's equations. 

In this work we extend some of these results to higher dimensional spacetimes. There are 
several reasons why one should now try to do it. It seems impossible to formulate in four 
dimensions a consistent theory which unifies gravity with the other forces in nature. Thus, 
most efforts in this direction have considered a higher dimensional arena for our universe, 
one example being string theories which have recently made some remarkable achievements. 
Moreover, recent investigations [9] propose the existence of extra dimensions in our Universe 
in order to solve the hierarchy problem, i.e., the huge difference between the electroweak 
and the Planck scale, mEw/^Pi ~ 10^^^. The fields of standard model would inhabit a 
4-dimensional sub-manifold, the brane, whereas the gravitational degrees of freedom would 
propagate throughout all dimensions. One of the most spectacular consequences of this 
scenario would be the production of black holes at the Large Hadron CoUider at CERN [10] 
(for recent relevant work related to this topic we refer the reader to [11-13]). Now, one of 
the experimental signatures of black hole production will be a missing energy, perhaps a 
large fraction of the center of mass energy [12]. This will happen because when the partons 
collide to form a black hole, some of the initial energy will be converted to gravitational 
waves, and due to the small amplitudes involved, there is no gravitational wave detector 
capable of detecting them, so they will appear as missing. Thus, the coUider could in fact 
indirectly serve as a gravitational wave detector. This calls for the calculation of the energy 
given away as gravitational waves when two high energy particles collide to form a black 
hole, which lives in all the dimensions. The work done so far on this subject [14, 15] in 
higher dimensions, is mostly geometric, and generalizes a construction by Penrose to find 
trapped surfaces on the union of two shock waves, describing boosted Schwarzschild black 
holes. On the other hand, there are clues [12, 13] indicating that a formalism described 
by Weinberg [16] to compute the gravitational energy radiated in the collision of two point 
particles, gives results correct to a order of magnitude when applied to the collision of two 
black holes. This formahsm assumes a hard collision, i.e., a coUision lasting zero seconds. 
It would be important to apply this formalism in higher dimensions, trying to see if there 
is agreement between both results. This will be one of the topics discussed in this paper. 
The other topic we study in this paper is the quadrupole formula in higher dimensions. Due 
to the difficulties in handling the wave tails in odd dimensions we concentrate our study in 
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even dimensions. 

This paper is organized as follows: In section II we linearize Einstein's equations in a flat 
£)-dimensional background and arrive at an inhomogeneous wave equation for the metric 
perturbations. The source free equations are analyzed in terms of plane waves, and then 
the general solution to the homogeneous equation is deduced in terms of the D-dimensional 
retarded Green's function. In section III we compute the D-dimensional quadrupole formula 
(assuming slowly moving sources), expressing the metric and the radiated energy in terms of 
the time-time component of the energy-momentum tensor. We then apply the quadrupole 
formula to two cases: a particle in circular motion in a generic background, and a particle 
falling into a D-dimensional Schwarzschild black hole. In section IV we consider the hard 
collision between two particles, i.e., the collision takes zero seconds, and introduce a cutoff 
frequency necessary to have meaningful results. We then apply to the case where one of 
the colliding particles is a black hole. We propose that this cutoff should be related to 
the gravitational quasinormal frequency of the black hole, and compute some values of the 
scalar quasinormal frequencies for higher dimensional Schwarzschild black holes, expecting 
that the gravitational quasinormal frequencies will behave in the same manner. Finally, we 
apply this formalism to compute the generation of gravitational radiation during black hole 
pair creation in four and higher dimensions, a result that has never been worked out, even 
for D — A. In our presentation we shall mostly follow Weinberg's [16] exposition. 

II. LINEARIZED Z)-DIMENSIONAL EINSTEIN'S EQUATIONS 

Due to the non-linearity of Einstein's equations, the treatment of the gravitational radia- 
tion problem is not an easy one since the energy-momentum tensor of the gravitational wave 
contributes to its own gravitational field. To overcome this difficulty it is a standard proce- 
dure to work only with the weak radiative solution, in the sense that the energy-momentum 
content of the gravitational wave is small enough in order to allow us to neglect its contri- 
bution to its own propagation. This approach is justified in practice since we expect the 
detected gravitational radiation to be of low intensity. 
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A. The inhomogeneous wave equation 



We begin this subsection by introducing the general background formahsm (whose details 
can be found, e.g., in [16]) that will be needed in later sections. Then we obtain the linearized 
inhomogeneous wave equation. 

Greek indices vary as 0,1, ■ ■ ■ , D — 1 and latin indices as 1, ■ ■ ■ , D—1 and our units are such 
that c = 1. Wc work on a D-dimensional spacctime described by a metric that approaches 
asymptotically the D-dimensional Minkowski metric 77^,^ = diag(— 1, +1, • • • , +1), and thus 
we can write 

fl'/.i/ = ^M^^ + V = 0, 1, • • • , i:) - 1 , (1) 

where hfj_„ is small, i.e., \hfj_^\ << 1, so that it represents small corrections to the fiat 
background. The exact Einstein field equations, G/^i, — SttQT^i, (with Q being the usual 
Newton constant), can then be written as 

with 

T'^^=r?'^V^7;^ + ^a/3)- (3) 

Here R^^\^ is the part of the Ricci tensor hnear in /i^^, R^^^'^a — rf^R^^^i^^, and t^v is the ef- 
fective energy-momentum tensor, containing contributions from T^,y, the energy-momentum 
tensor of the matter source, and t^,^ which represents the gravitational contribution. The 
pseudo-tensor tfjtv contains the difference between the exact Ricci terms and the Ricci terms 
linear in h^^, 

RfJ-i^ ~ 7,9lJ.fR a ~ R^ \v + TzVlJ-i^R^ ^ a 



The Bianchi identities imply that r^^ is locally conserved. 



(4) 



d^T'^'' = . (5) 

Introducing the cartesian coordinates = (i,x) with x = x\ and considering a D — 1 
volume V with a boundary spacelike surface S with dimension D — 2 whose unit exterior 
normal is n, eq. (5) yields 



d 
~dt 



[ d^-^x T^'^ = - / d^-^x mT'" . (6) 
Jv Js 
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This means that one may interpret 

p''= [ d^-ixT°" (7) 
Jv 

as the total energy-momentum (pseudo)vector of the system, including matter and gravita- 
tion, and r**^ as the corresponding flux. Since the matter contribution is contained in f^" , 
the flux of gravitational radiation is 



Flux= / (i^-^xn^r. (8) 

5 



In this context of linearized general relativity, we neglect terms of order higher than 
the first in h^i, and all the indices are raised and lowered using 77^*^. We also neglect the 
contribution of the gravitational energy-momentum tensor t^,^ (i.e., << \T^u\) since 
from (4) we see that tf^i, is of higher order in h^i,. Then, the conservation equations (5) yield 

d^T"" = . (9) 

In this setting and choosing the convenient coordinate system that obeys the harmonic (also 
called Lorentz) gauge conditions, 

2a^/i\ = a,r, (10) 

(where — d/dx^), the first order Einstein field equations (2) yield 

□V = -167r6?^^. , (11) 

where □ = ri^^d^di, is the D-dimensional Laplacian, and S^i, will be called the modified 
energy-momentum tensor of the matter source. Eqs. (11) and (12) subject to (10) allow us 
to find the gravitational radiation produced by a matter source Sfj^v 



B. The plane wave solutions 

In vacuum, the linearized equations for the gravitational field are R^^\v — or, equiva- 
lently, the homogeneous equations \^hfj^v — 0, subjected to the harmonic gauge conditions 
(10). The solutions of these equations, the plane wave solutions, are important since the 
general solutions of the inhomogeneous equations (10) and (11) approach the plane wave 
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solutions at large distances from the source. Setting ka — (— k) with cu and k being 
respectively the frequency and wave vector, the plane wave solutions can be written as a 
linear superposition of solutions of the kind 

K.{t, x) = e,, e'^"^" + el, e-'^-^" , (13) 

where e^i, — 61,^ is called the polarization tensor and * means the complex conjugate. These 
solutions satisfy eq. (11) with S/j^i, = if kak'^ — 0, and obey the harmonic gauge conditions 
(10) if 2A;^e^, = A;,e'^^. 

An important issue that must be addressed is the number of different polarizations that 
a gravitational wave in D dimensions can have. The polarization tensor e^j,, being symmet- 
ric, has in general D{D + l)/2 independent components. However, these components are 
subjected to the D harmonic gauge conditions that reduce by D the number of independent 
components. In addition, under the infinitesimal change of coordinates x''^ — x'^ + ^'^{x), the 
polarization tensor transforms into e^^, — e^i, — 81,$,^ — d/^^^- Now, e^, and e^i^ describe the 
same physical system for arbitrary values of the D parameters $,'^{x). Therefore, the number 
of independent components of e^,y, i.e., the number of polarization states of a gravitational 
wave in D dimensions is D{D + l)/2 — D — D = D[D — 3)/2. From this computation we 
can also see that gravitational waves are present only when D > 2>. Therefore, from now on 
we assume D > 3 whenever we refer to D. In what concerns the helicity of the gravitational 
waves, for arbitrary D the gravitons are always spin 2 particles. 

To end this subsection on gravitational plane wave solutions, we present the average 
gravitational energy-momentum tensor of a plane wave, a quantity that will be needed 
later. Notice that in vacuum, since the matter contribution is zero (T^j, = 0), we cannot 
neglect the contribution of the gravitational energy-momentum tensor t^,. From eq. (4), 
and neglecting terms of order higher than h?, the gravitational energy- momentum tensor of 
a plane wave is given by 

1 



(14) 



and through a straightforward calculation (see e.g. [16] for details) we get the average 
gravitational energy-momentum tensor of a plane wave. 



af} * |2 



(15) 
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C. The D-dimensional retarded Green's function 



The general solution to the inhomogeneous differential equation (11) may be found in the 
usual way in terms of a Green's function as 

hni,{t,x.) — —IGnQ J dt' J d^~^x!Si^i,{t\x')G{t — t',x — x')+hom.ogeneoussohjLtions, (16) 

where the Green's function G{t — t', x — x') satisfies 

rj'^'d^.d^^Git - i', X - x') = S{t - t')5{yi - x') , (17) 

where 5{z) is the Dirac delta function. In the momentum representation this reads 

G{t, x) = -— -- / d'>-'ke^^ - / du^^— , (18) 



(27r)^ J J uj'^-k'^' 

where fc^ = /cf + /c| + ... + To evaluate this, it is convenient to perform the fc-integral by 

using spherical coordinates in the (D — l)-dimensional /c-space. The required transformation, 
along with some useful formulas which shall be used later on, is given in Appendix A. The 
result for the retarded Green's function in those spherical coordinates is 

G'''^ X) = -j^^^, X ;:(^ / /c(^-^)/Vp_3)/2(/cr) Mkt)dk , (19) 

where r'^ — %{ + x\-\- ... + and 6(t) is the Heaviside function defined as 

. X f 1 if ^ > 

m = { (20) 

[ if i < . 

The function J(d-3)/2(^'") is a Bcsscl function [17, 18]. The structure of the retarded Green's 
function will depend on the parity of as we shall see. This dependence on the parity, 
which implies major differences between even and odd spacetime dimensions, is connected 
to the structure of the Bessel function. For even D, the index of the Bessel function is 
semi-integer and then the Bessel function is expressible in terms of elementary functions, 
while for odd D this does not happen. A concise explanation of the difference between 
retarded Green's function in even and odd and the physical consequences that entails is 
presented in [19] (see also [20-22]). A complete derivation of the Green's function in higher 
dimensional spaces may be found in Hassani [23]. The result is 



G"*(i,x) = ± 



d ' 


(D-4)/2 


'6{t-ry 


2Trrdr 




r 



D even. (21) 
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d 

2'Krdr 



D odd. (22) 



It is sometimes convenient to work with the Fourier transform (in the time coordinate) of 
the Green's function. One finds [23] an analytical result independent of the parity of D 

X) = ^(^^^y(^ (-) i^?.-3)/2(-) , (23) 

where Hl{z) is a modified Bessel function [17, 18]. Of course, the different structure of the 
Green's function for different D is again embodied in these Bessel functions. Equations (21) 
and (23), are one of the most important results we shall use in this paper. For D = A (21) 
obviously reproduce well known results [23]. Now, one sees from eq. (21) that although there 
are delta function derivatives on the even-D Green's function, the localization of the Green's 
function on the light cone is preserved. However, eq. (22) tells us that the retarded Green's 
function for odd dimensions is non-zero inside the light cone. The consequence, as has been 
emphasized by different authors [19, 22, 24], is that for odd D the Huygens principle does 
not hold: the fact that the retarded Green's function support extends to the interior of the 
light cone implies the appearance of radiative tails in (16). In other words, we still have a 
propagation phenomenum for the wave equation in odd dimensional spacetimes, in so far as 
a localized initial state requires a certain time to reach a point in space. Huygens principle 
no longer holds, because the effect of the initial state is not sharply limited in time: once 
the signal has reached a point in space, it persists there indefinitely as a reverberation. 

This fact coupled to the analytic structure of the Green's function in odd dimensions make 
it hard to get a grip on radiation generation in odd dimensional spacetimes. Therefore, from 
now on we shall focus on even dimensions, for which the retarded Green's function is given 
byeq. (21). 

D. The even D-dimensional retarded solution in the wave zone 

The retarded solution for the metric perturbation h^^, obtained by using the retarded 
Green's function (21) and discarding the homogeneous solution in (16) will be given by 

x) = 167r^? j dt' J d'^-^x'S^^it', ^')G''\t - i', X - x') , (24) 

with C*^*(t — t', X — x') as in eq. (21). For D = 4 for example one has 

G-%.) = 1^-^, D^A, (25) 
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which is the well known result. For D = 6, we have 

where the S'{t — r) means derivative of the Dirac delta function with respect to its argument. 
For D — 8, we have 

We see that in general even-D dimensions the Green's function consists of inverse integer 
powers in r, spanning all values between ^(d-2)/2 and ^7=3, including these ones. Now, the 
retarded solution is given by eq. (24) as a product of the Green's function times the modified 
energy-momentum tensor Sfj^^- The net result of having derivatives on the delta functions 
is to transfer these derivatives to the energy-momentum tensor as time derivatives (this can 
be seen by integrating (24) by parts in the ^-integral). 

A close inspection then shows that the retarded field possesses a kind of peeling property 
in that it consists of terms with different fall off at infinity. Explicitly, this means that the 
retarded field will consist of a sum of terms possessing all integer inverse powers in r between 
and D — 3. The term that dies off more quickly at infinity is the typically a static 

term, since it comes from the Laplacian. As a matter of fact this term was already observed 
in the higher dimensional black hole by Tangherhni [25] (see also Myers and Perry [26]). 
We will see that the term faUing more slowly, the one that goes like ^(0-2) /2 : gives rise to 
gravitational radiation. It is well defined, in the sense that the power crossing sufficiently 
large hyperspheres with different radius is the same, because the volume element goes as 
r^"'^ and the energy as ~ 

In radiation problems, one is interested in finding out the field at large distances from 
the source, r >> A, where A is the wavelength of the waves, and also much larger than the 
source's dimensions R. This is defined as the wave zone. In the wave zone, one may neglect 
all terms in the Green's function that decay faster than (d-2)/2 ■ So, in the wave zone, we 
find 



1 (2^) 



J rf^-V5^,(t- |x-x'|,x') 



(28) 



(27rr)(^-2)/2 * 

where ^ stands for the '^^th derivative with respect to time. For D = 4 eq. (28) yields 
the standard result [16]: 

V(^,x) = -^y"c^^-V^^,(i-|x-x'|,x'), D^4. (29) 
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To find the Fourier transform of the metric, one uses the representation (23) for the Green's 
function. Now, in the wave zone, the Green's function may be simphfied using the asymptotic 
expansion for the Bessel function [18] 



i^f.-3)/.(-r) ~ ^ ^^4^^^-^ , a;r ^ oo. (30) 

This yields 

x) = -^^-^^^a;(^-^)/^e- J d^-^^'S,^{u, x') . (31) 

This could also have been arrived at directly from (28), using the rule time derivative — > —icu 
for Fourier transforms. Equations (28) and (31) are one of the most important results derived 
in this paper, and will be the basis for all the subsequent section. Similar equations, but not 
as general as the ones presented here, were given by Chen, Li and Lin [27] in the context of 
gravitational radiation by a rolling tachyon. 

To get the energy spectrum, we use (12) yielding 

III. THE EVEN Z)-DIMENSIONAL QUADRUPOLE FORMULA 

A. Derivation of the even D-dimensional quadrupole formula 

When the velocities of the sources that generate the gravitational waves are small, it 
is sufficient to know the component of the gravitational energy-momentum tensor in 
order to have a good estimate of the energy they radiate. In this subsection, we will deduce 
the D-dimensional quadrupole formula and in the next subsection we will apply it to (1) a 
particle in circular orbit and (2) a particle in free fall into a D-dimensional Schwarzschild 
black hole. 

We start by recalling that the Fourier transform of the energy-momentum tensor is 

T^^iu, k) = y d^-'x'e-'^ ''' J dt e^' T^^it, x) + c.c. , (33) 

where c.c. means the complex conjugate of the preceding term. Then, the conservation 
equations (9) for T^''{t,x) applied to eq. (33) yield k^T^i,{u;,'k) — 0. Using this last result 
we obtain Too{uj, k) = Tji{uj, k) and Toi{uj, k) = —k^ Tji{uj, k), where k = k/o;. We can 
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then write the energy spectrum, eq. (32), as a function only of the spacehke components of 
r'^'^(u;,k), 

'^'^ 2g-^^K^j^{k)T*'^{uj,k)T^^{u,k) , (34) 



dudn (27r)^-4 
where 

■^ij,lm{k) = SiiSjra ~ '^kjkmSu + ^ {^^ijSim + kikm^ij + kikjSim^ + — ~kikjkikm ■ (35) 

At this point, we make a new approximation (in addition to the wave zone approximation) 
and assume that uiR << 1, where R is the source's radius. In other words, we assume that 
the internal velocities of the sources are small and thus the source's radius is much smaller 
than the characteristic wavelength ~ l/cu of the emitted gravitational waves. Within this 
approximation, one can set e"'"*'^' ~ 1 in eq. (33) (since R — Ix'jmax)- Moreover, after 
a straightforward calculation, one can also set in eq. (34) the approximation T'-'(a;,k) ~ 
— {LJ^/2)Dij{ui), where 

Dij{u) = y"(i^-'xx'x^T°°(c^,x). (36) 



Finally, using 



/ 



/ 



dfl£i-2kikj — — -Sij , 



3Q^ 2 

dVtD-2kikjkikm = £)2 ~i i^ij^lrn + 5ii5jm + 5im5jl) , (37) 



where f2D_2 is the {D — 2)-dimensional solid angle defined in (A5), we obtain the D- 
dimensional quadrupole formula 



dE _ 22-^7r-(-°-5)/2^ {D - 3)D 



D+2 



(D-l)m(a;)A,(a;)-|A^(^)|' , (38) 



dw r[{D - 1) /2] (L>2 -1){D- 2) 

where the Gamma function T[z] is defined in Appendix A. As the dimension D grows it is 
seen that the rate of gravitational energy radiated increases as u^^"^. Sometimes it will be 
more useful to have the time rate of emitted energy 

(39) 



dt r[{D-l)/2]{D'^-l){D-2) 



For D — 4, eq. (39) yields the well known result [16] 



dE_g 
dt ~ 5 



d!D;^{mD,,{t)-\\c^Du{t)\']. (40) 
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B. Applications of the quadrupole formula: test particles in a background geom- 
etry 

The quadrupole formula has been used successfully in almost all kind of problems in- 
volving gravitational wave generation. By successful we mean that it agrees with other 
more accurate methods. Its simplicity and the fact that it gives results correct to within a 
few percent, makes it an invaluable tool in estimating gravitational radiation emission. We 
shall in the following present two important examples of the application of the quadrupole 
formula. 



1. A particle in circular orbit 

The radiation generated by particles in circular motion was perhaps the first situation to 
be considered in the analysis of gravitational wave generation. For orbits with low frequency, 
the quadrupole formula yields excellent results. As expected it is difficult to find in nature a 
system with perfect circular orbits, they will in general be elliptic. In this case the agreement 
is also remarkable, and one finds that the quadrupole formalism can account with precision 
for the increase in period of the pulsar PSR 1913+16, due to gravitational wave emission [28]. 
In four dimensions the full treatment of elliptic orbital motion is discussed by Peters [29]. In 
dimensions higher than four, it has been shown [25] that there are no stable geodesic circular 
orbits, and so geodesic circular motion is not as interesting for higher D. For this reason, and 
also because we only want to put in evidence the differences that arise in gravitational wave 
emission as one varies the spacetime dimension D, we will just analyze the simple circular, 
not necessarily geodesic motion, to see whether the results are non-trivially changed as one 
increases D. Consider then two bodies of equal mass m in circular orbits a distance / apart. 
Suppose they revolve around the center of mass, which is at //2 from both masses, and that 
they orbit with frequency lu in the x — y plane. A simple calculation [29, 30] yields 

Dxx = cos{2ujt) + const , (41) 

Dyy — Dxx 1 ('^2) 

Dxy = —7- sin(2a;i) + const , (43) 
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independently of the dimension in which they are imbedded and with all other components 
being zero. We therefore get from eq. (39) 



^ ^ ^2^4^Z.+2 .44. 

dt 7r(-D-5)/2r[(L> - 1) /2] {D + 1){D- 2) " ^ ^ 

For D — A one gets 

-m'Z^o;^ , (45) 



rft 5 

which agrees with known results [29, 30]. Eq. (44) is telling us that as one climbs up in 
dimension number D, the frequency effects gets more pronounced. 

2. A particle falling radially into a higher dimensional Schwarzschild black hole 

As yet another example of the use of the quadrupole formula eq. (39) we now calculate 
the energy given away as gravitational waves when a point particle, with mass m falls into 
a D-dimensional Schwarzschild black hole, a metric first given in [25]. Historically, the 
case of a particle falling into a D — 4 Schwarzschild black hole was one of the first to 
be studied [31, 32] in connection with gravitational wave generation, and later served as a 
model calculation when one wanted to evolve Einstein's equations fully numerically [33, 34] . 
This process was first studied [32] by solving numerically ZeriUi's [31] wave equation for a 
particle at rest at infinity and then falling into a Schwarzschild black hole. Davis et al [32] 
found numerically that the amount of energy radiated to infinity as gravitational waves was 

2 

Ai?num = 0.01^, where m is the mass of the particle falling in and M is the mass of the 
black hole. It is found that the D — 4 quadrupole formula yields [35] AE^^^ — 0.019^, so 
it is of the order of magnitude as that given by fully relativistic numerical results. Despite 
the fact that the quadrupole formula fails somewhere near the black hole (the motion is not 
slow, and the background is certainly not fiat), it looks like one can get an idea of how much 
radiation will be released with the help of this formula. Based on this good agreement, we 
shall now consider this process but for higher dimensional spacetimes. The metric for the 
D-dimensional Schwarzschild black hole in {t, r, 9i, 62, .., ^0-2) coordinates (see Appendix A) 
is 

(46) 
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Consider a particle falling along a radial geodesic, and at rest at infinity. Then, the geodesic 
equations give 

dr IGnQM 1 

where we make the flat space approximation t — t. We then have, in these coordinates, 
Dii — r^, and all other components vanish. Prom (39) we get the energy radiated per 
second, which yields 

lit - r[(D-i)/2](D2-i) ^1^* ^^^1 ' 

We can perform the derivatives and integrate to get the total energy radiated. There is a 
slight problem though, where do we stop the integration? The expression for the energy 
diverges at r = but this is no problem, as we know that as the particle approaches the 
horizon, the radiation will be infinitely red-shifted. Moreover, the standard picture [35] is 
that of a particle falling in, and in the last stages being frozen near the horizon. With this 
in mind we integrate from r = oo to some point near the horizon, say r = 6 x r_|_, where r+ 
is the horizon radius and b is some number larger than unit, and we get 

.D(D — 2)n /Q n2wo 



(50) 



where 

(3 - L»)2(5 - Dfiy - 3Dy{8 - 4D)2(9 - hD)\..{D/2 + 4 - /2f 
~ T[{D - l)/2f{D -!)(£'+ 1){D + 3) 

To understand the effect of both the dimension number D and the parameter h on the total 
energy radiated according to the quadrupole formula, we list in Table 1 some values A£' for 
different dimensions, and h between 1 and 1.3. 

The parameter h is in fact a measure of our ignorance of what goes on near the black 
hole horizon, so if the energy radiated doesn't vary much with h it means that our lack of 
knowledge doesn't affect the results very much. For D — A that happens indeed. Putting 
6=1 gives only an energy 2.6 times larger than with h = 1.3, and still very close to the fully 
relativistic numerical result of 0.01^. However as we increase the effect of h increases 
dramatically. For D = 12 for example, we can see that a change in h from 1 to 1.3 gives 
a corresponding change in l\E of 3 x 10^ to 0.0665. This is 8 orders of magnitude lower! 
Since there is as yet no Regge- Wheeler- Zerilli [31, 36] wavefunction for higher dimensional 
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TABLE I: The energy 


radiated by 


a particle falling from rest into a higher 


dimensional 


Schwarzschild black hole, 


as a function of dimension. The integration is stopped at b 


X r+ where 


r+ is the horizon radius. 










D 


6=1: 


b= 1.2: 


b = 1.3: 


4 


0.019 


0.01 


0.0076 


6 


0.576 


0.05 


0.0167 


8 


180 


1.19 


0.13 


10 


24567 


6.13 


0.16 


12 


3.3 X 10^ 


14.77 


0.0665 



Schwarzschild black holes, there are no fully relativistic numerical results to compare our 
results with. Thus D = 4 is just the perfect dimension to predict, through the quadrupole 
formula, the gravitational energy coming from collisions between particles and black holes, 
or between small and massive black holes. It is not a problem related to the quadrupole 
formalism, but rather one related to -D. A small change in parameters translates itself, 
for high D, in a large variation in the final result. Thus, as the dimension D grows, the 
knowledge of the cutoff radius 6 x r+ becomes essential to compute accurately the energy 
released. 



IV. INSTANTANEOUS COLLISIONS IN EVEN D-DIMENSIONS 

In general, whenever two bodies collide or scatter there will be gravitational energy re- 
leased due to the changes in momentum involved in the process. If the collision is hard 
meaning that the incoming and outgoing trajectories have constant velocities, there is a 
method first envisaged by Weinberg [16, 37], later explored in [38] by Smarr to compute 
exactly the metric perturbation and energy released. The method is valid for arbitrary ve- 
locities (one will still be working in the linear approximation, so energies have to be low). 
Basically, it assumes a collision lasting for zero seconds. It was found that in this case the 
resulting spectra were fiat, precisely what one would expect based on one's experience with 
electromagnetism [39], and so to give a meaning to the total energy, a cutoff frequency is 
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needed. This cutoff frequency depends upon some physical cutoff in the particular problem. 
We shall now generalize this construction for arbitrary dimensions. 



A. Derivation of the Radiation Formula in Terms of a Cutoff for a Head-on Col- 
lision 

Consider therefore a system of freely moving particles with L>-momenta Pj'", energies Ei 
and l)-velocities v, which due to the coUision change abruptly ait — 0, to corresponding 
primed quantities. For such a system, the energy-momentum tensor is 

T^'^it, v) = ^ ^'^''"'(^ - ^tM-t) + -^s^'-'i^' - Vt)e{t) , (51) 

from which, using cqs. (31) and (32) one can get the quantities /i^j^ and also the radiation 
emitted. Let us consider the particular case in which one has a head-on collision of two 
particles, particle 1 with mass mi and Lorentz factor 71, and particle 2 with mass m — m2 
with Lorentz factor 72, colliding to form a particle at rest. Without loss of generahty, one 
may orient the axis so that the motion is in the (xd_i,xd) plane, and the xd axis is the 
radiation direction (see Appendix A). We then have 

Pi =7imi(l,0,0, ...,'f;isin^i,'t;iCOS^i) ; P^' = (£;i, 0, 0, 0, 0) (52) 
P2 = 72?«2(1,0,0,...,-V2sin^i,-V2cos^i) ; P^ = (E^, 0, 0, 0, 0). (53) 

Momentum conservation leads to the additional relation ^\m\Vi — 72m2t'2. Replacing (52) 
and (53) in the energy-momentum tensor (51) and using (32) we find 

(PE ^ 2g D-3 -ffmlvj{vi + V2fsm9,^ ^-4 /..n 

dujdQ (27^)^-2p>-2(l-^;lCOS^l)2(l + ^;2COS^l)2 ' ^ ^ 

We see that the for arbitrary (even) D the spectrum is not fiat. Flatness happens only for 
D — A. For any D the total energy, integrated over all frequencies would diverge so one 
needs a cutoff frequency which shall depend on the particular problem under consideration. 
Integrating (54) from = to the cutoff frequency Uc we have 

dE ^ 2g 1 -ffmjvf{v^ + V2fsme,* 

dn (27r)^-2p>-2(l-^;iCOS^i)2(l + ^;2COS^i)2 ^ '^'^ ' ^ ^ 

Two limiting cases are of interest here, namely (i) the collision between identical particles 
and (ii) the collision between a light particle and a very massive one. In case (i) replacing 
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mi = 777,2 — m, vi — V2 — V, eq. (55) gives 

(27r)^-2L»-2(l-v2cos2^i)2 ^'^'^ ■ ^ ^ 

In case (ii) considering mi7i = 7777 << 777272, f 1 = f >> f2, eq. (55) yields 

dE _ 2g 1 7Wt;Sing/ ^,3 

~ (27r)^-2L>-2(l-^;cos^i)2 ^ ' ^ ^ 

Notice that the technique just described is expected to break down if the velocities involved 
are very low, since then the collision would not be instantaneous. In fact a condition for 
this method to work would can be stated 

Indeed, one can see from eq. (55) that if f ^ 0, ^ — 0, even though wc know (see 
Subsection (IV A)) that AE ^ 0. In any case, if the velocities are small one can use the 
quadrupole formula instead. 



B. Applications: The Cutoff Frequency when one of the Particles is a Black Hole 
and Radiation from Black Hole Pair Creation 

1. The Cutoff Frequency when one of the Head-on Colliding Particles is a Black Hole 

We shall now restrict ourselves to the case (ii) of last subsection, in which at least one 
of the particles participating in the collision is a massive black hole, with mass M » m 
(where we have put r?7-i = r?7 and r?72 = M). Formulas (55)- (57) are useless unless one is able 
to determine the cutoff frequency ujc present in the particular problem under consideration. 
In the situation where one has a small particle colliding at high velocities with a black 
hole, it has been suggested by Smarr [38] that the cutoff frequency should be Uc ~ 1/2M, 
presumably because the characteristic collision time is dictated by the large black hole whose 
radius is 2M. Using this cutoff he finds 

2 

777 

A£;s„,arr ~ 0.27^ — . (58) 

The exact result, using a relativistic perturbation approach which reduces to the numerical 
integration of a second order differential equation (the Zerilli wavefunction) , has been given 
by Cardoso and Lemos [13], as 

2 

777 

A£;exa.t = 0.2672— . (59) 
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m-y - dco 



0.26 



0.075 
0.03 




FIG. 1: The energy spectra as a function of the angular number for a highly relativistic particle 
falling into a D = 4 Schwarzschild black hole [13]. The particle begins to fall with a Lorentz factor 
7. Notice that for each / there is a cutoff frequency coci which is equal to the quasinormal frequency 
a;QN after which the spectrum decays exponentially. So it is clearly seen that a;QN works as a cutoff 
frequency. The total energy radiated is a given by a sum over I, which is the same as saying that 
the effective cutoff frequency is given by a weighted average of the various . 

This is equivalent to saying that Uc = ~ TgIm' looks like the cutoff is 

indeed the inverse of the horizon radius. However, in the numerical work by Cardoso and 
Lemos, it was found that it was not the presence of an horizon that contributed to this 
cutoff, but the presence of a potential barrier V outside the horizon. By decomposing 
the field in tensorial spherical harmonics with index I standing for the angular quantum 
number, we found that for each I, the spectrum is indeed fiat (as predicted by eq. (55) for 
D = 4:), until a cutoff frequency ujc^ which was numerically equal to the lowest gravitational 
quasinormal frequency cjqn. For u > uj^ the spectrum decays exponentially. This behavior is 
illustrated in Fig. 1. The quasinormal frequencies [40] are those frequencies that correspond 
to only outgoing waves at infinity and only ingoing waves near the horizon. As such the 
gravitational quasinormal frequencies will in general have a real and an imaginary part, 
the latter denoting gravitational wave emission and therefore a decay in the perturbation. 
There have been a wealth of works dwelling on quasinormal modes on asymptotically fiat 
spacetimes [40], due to its close connection with gravitational wave emission, and also on 
non-asymptotically flat spacetimes, like asymptotically anti-de Sitter [41] or asymptotically 
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de Sitter [42] spacetimes, mainly due to the AdS/CFT and dS/CFT [43] correspondence 
conjecture. We argue here that it is indeed the quasinormal frequency that dictates the 
cutoff, and not the horizon radius. For = 4 it so happens that the weighted average of 
uOf.^ is ^^"Ip, which, as we said, is quite similar to r_|_ = j^. The reason for the cutoff being 
dictated by the quasinormal frequency can be understood using some WKB intuition. The 
presence of a potential barrier outside the horizon means that waves with some frequencies 
get reflected back on the barrier while others can cross. Frequencies such that o;^ is lower 
than the maximum barrier height V^ax will be reflected back to infinity where they will be 
detected. However, frequencies a;^ larger than the maximum barrier height cross the barrier 
and enter the black hole, thereby being absorbed and not contributing to the energy detected 
at infinity. So only frequencies cu^ lower than this maximum barrier height are detected at 
infinity. It has been shown [44] that the gravitational quasinormal frequencies are to first 
order equal to the square root of the maximum barrier height. In view of this picture, and 
considering the physical meaning of the cutoff frequency, it seems quite natural to say that 
the cutoff frequency is equal to the quasinormal frequency. If the frequencies are higher than 
the barrier height, they don't get reflected back to infinity. This discussion is very important 
to understand how the total energy varies with the number D of dimensions. In fact, if we 
set a;c ~ V^i ^^^^ total energy radiated decreases rapidly with the dimension 

number, because r+ increases rapidly with the dimension. This conflicts with recent results 
[14, 15], which using shock waves that describe boosted Schwarzschild black holes, and 
searching for apparent horizons, indicate an increase with D. So, we need the gravitational 
quasinormal frequencies for higher dimensional Schwarzschild black holes. To arrive at an 
wave equation for gravitational perturbations of higher dimensional Schwarzschild black 
holes, and therefore to compute its gravitational quasinormal frequencies, one needs to 
decompose Einstein's equations in D-dimensional tensorial harmonics, which would lead to 
some quite complex expressions. It is not necessary to go that far though, because one 
can get an idea of how the gravitational quasinormal frequencies vary by searching for the 
quasinormal frequencies of scalar perturbations, and scalar quasinormal frequencies are a 
lot easier to find. One hopes that the scalar frequencies will behave with D in the same 
manner as do the gravitational ones. Scalar perturbations in D-dimensional Schwarzschild 
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TABLE II: The lowest scalar quasinormal frequencies for spherically symmetric (/ = 0) scalar 
perturbations of higher dimensional Schwarzschild black holes, obtained using a WKB method 
[44] . Notice that the real part of the quasinormal frequency is always the same order of magnitude 
as the square root of the maximum barrier height. We show also the maximum barrier height as 
well as the horizon radius as a function of dimension D. The mass M of the black hole has been 
set to 1. 



D 


RefwQjv]: 


ImfwQjv] : 


\/Vmax: 


l/r+: 


4 


0.10 


-0.12 


0.16 


0.5 


6 


1.033 


-0.713 


1.441 


1.28 


8 


1.969 


-1.023 


2.637 


1.32 


10 


2.779 


-1.158 


3.64 


1.25 


12 


3.49 


-1.202 


4.503 


1.17 



spacetimes obey the wave equation (consult [45] for details) 

d'^(p{uj, r 



drl 

The potential V{r) appearing in equation (60) is given by 
V{r) = f{r) 



+ -V{r)](l){u,r) ^0. (60) 



a_ ^ {D-2){D-4)f{r) ^ {D - 2)/'(r) 



(61) 



where a — l{l+D—S) is the eigenvalue of the Laplacian on the hypersphere S^'"^, the tortoise 
coordinate n is defined as |^ = /(r) = (l - ^^i||M__^j, and f'{r) = We have 

found the quasinormal frequencies of spherically symmetric (/ = 0) scalar perturbations, by 
using a WKB approach developed by Schutz, Will and collaborators [44]. The results are 
presented in Table 2, where we also show the maximum barrier height of the potential in 
eq. (61), as well as the horizon radius. 

The first thing worth noticing is that the real part of the scalar quasinormal frequency 



is to first order reasonably close to the square root of the maximum barrier height VKnax, 
supporting the previous discussion. Furthermore, the scalar quasinormal frequency grows 
more rapidly than the inverse of the horizon radius ^ as one increases D. In fact, the scalar 
quasinormal frequency grows with D while the horizon radius r+ gets smaller. Note that 
from pure dimensional arguments, for fixed D, a; oc The statement here is that the 
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constant of proportionality depends on the dimension more explicitly it grows with 
and can be found from Table 2. Assuming that the gravitational quasinormal frequencies 
will have the same behavior (and some very recent studies [46] relating black hole entropy 

and damped quasinormal frequencies seem to point that way), the total energy radiated will 
during high-energy collisions does indeed increase with as some studies [14, 15] seem to 
indicate. 

2. The gravitational energy radiated during black hole pair creation 

As a new application of this instantaneous collision formalism, we will now consider the 
gravitational energy released during the quantum creation of pairs of black holes, a process 
which as far as we know has not been analyzed in the context of gravitational wave emission, 
even for D = 4. It is well known that vacuum quantum fluctuations produce virtual electron- 
positron pairs. These pairs can become real [47] if they are pulled apart by an external 
electric field, in which case the energy for the pair materialization and acceleration comes 
from the external electric field energy. Likewise, a black hole pair can be created in the 
presence of an external field whenever the energy pumped into the system is enough in 
order to make the pair of virtual black holes real (see Dias [48] for a review on black hole 
pair creation). If one tries to predict the spectrum of radiation coming from pair creation, 
one expects of course a spectrum characteristic of accelerated masses but one also expects 
that this follows some kind of signal indicating pair creation. In other words, the process of 
pair creation itself, which involves the sudden creation of particles, must imply emission of 
radiation. It is this phase we shall focus on, forgetting the subsequent emission of radiation 
caused by the acceleration. 

Pair creation is a pure quantum-mechanical process in nature, with no classical explana- 
tion. But given that the process does occur, one may ask about the spectrum and intensity 
of the radiation accompanying it. The sudden creation of pairs can be viewed for our pur- 
poses as an instantaneous creation of particles (i.e., the time reverse process of instantaneous 
collisions), the violent acceleration of particles initially at rest to some final velocity in a 
very short time, and the technique described at the beginning of this section applies. This 
is quite similar to another pure quantum-mechanical process, the beta decay. The electro- 
magnetic radiation emitted during beta decay has been computed classically by Chang and 
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Falkoff [49] and is also presented in Jackson [39]. The classical calculation is similar in all 
aspects to the one described in this section (the instantaneous collision formalism) assuming 
the sudden acceleration to energies E oi a charge initially at rest, and requires also a cutoff 
in the frequency, which has been assumed to be given by the uncertainty principle cj^ ~ • 
Assuming this cutoff one finds that the agreement between the classical calculation and the 
quantum calculation [49] is extremely good (specially in the low frequency regime), and 
more important, was verified experimentally. Summarizing, formula (56) also describes the 
gravitational energy radiated when two black holes, each with mass m and energy E form 
through quantum pair creation. The typical pair creation time can be estimated by the 
uncertainty principle Tcreation ~ ^^"^ ^^^^ cutoff frequency as 

■ ' (62) 



^creation ^ 

Here we would like to draw the reader's attention to the fact that the units of Planck's 
constant H change with dimension number D: according to our convention of setting c—1 

D-2 

the units of h are [M\D-i . With this cutoff, we find the spectrum of the gravitational 

radiation emitted during pair creation to be given by (54) with mi = m2 and Vi = V2 (we 
are considering the pair creation of two identical black holes) : 

(fE Sg D-3 7'^^iVsin^^/ 0-4 /«o^ 

X u! , (63) 



dudn (27r)^-2 L> - 2 (1 - cos^ 61)^ 

and the total frequency integrated energy per solid angle is 

dE _ Sg v^sin^i^ (7717)^-^ 

dn ~ (27r)^-2(D - 2) (1 - ^2 cos2 ^1)2 ^ h^-^ ■ 

For example, in four dimensions and for pairs with v ^ 1 one obtains 

dE 46? 2 2 ff,r^ 
— = — 7 m , 65) 

dcu TT 

and will have for the total energy radiated during production itself, using the cutoff frequency 
(62) 

AB=!?l!^. (66) 
TT n 

This could lead, under appropriate numbers of m and 7 to huge quantities. Although one 
cannot be sure as to the cutoff frequency, and therefore the total energy (66), it is extremely 
likely that, as was verified experimentally in beta decay, the zero frequency limit, eq. (65), 
is exact. 
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V. SUMMARY AND DISCUSSION 



We have developed the formahsm to compute gravitational wave generation in higher 
D dimensional spacetimes, with D even. Several examples have been worked out, and one 
cannot help the feeling that our apparently four dimensional world is the best one to make 
predictions about the intensity of gravitational waves in concrete situations, in the sense that 
a small variation of parameters leads in high D to a huge variation of the energy radiated. 
A lot more work is still needed if one wants to make precise predictions about gravitational 
wave generation in D dimensional spacetimes. For example, it would be important to find 
a way to treat gravitational perturbations of higher dimensional Schwarzschild black holes. 
One of the examples worked out, the gravitational radiation emitted during black hole pair 
creation, had not been previously considered in the literature, and it seems to be a good 
candidate, even in D = 4, to radiate intensely through gravitational waves. 
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APPENDIX A: SPHERICAL COORDINATES IN {D - l)-DIMENSIONS 



In this appendix we list some important formulas and results used throughout this pa- 
per. We shall first present the transformation mapping a (D — 1) cartesian coordinates, 
(xi, X2, X3, xd-i) onto {D — 1) spherical coordinates, (r, 9i, 62, 00-2) ■ The transforma- 
tion reads 

Xi — r sin 9i sin ^2--- sin^D_2 
X2—r sin ^1 sin ^2 ■ ■ ■ sin 6d-z cos 9d-2 



Xi — r sin 9i sin ^2 ■ ■ ■ sin ^ d-i-i cos 9 o-i 

xd-i = r cos 6*1 

The Jacobian of this transformation is 

J = sin 9^-^ sin 9^-\.. sin 9^-'-\.. sin 9d-3 , (Al) 
and the volume element becomes 

d^^^x = Jdrd9id92...d9D-2 = r^''^drdVtD-2 , (A2) 

where 

dnD-2 = sin 9^-^ sin 9^-^... sin 9D-zd9id92...d9D-2- (A3) 
is the element of the {D — \) dimensional solid angle. Finally, using 



!\„m- + l)/2] ,,,, 

y„ -^rl(n + 2)/2l' 

this yields 

2t^{D-1)/2 

^^-^ = r[p-i)/2]- ^"^^^ 

Here, V[z\ is the Gamma function, whose definition and properties are listed in [18]. In this 
work the main properties of the Gamma function which were used are V[z + 1] = ^^[A ^-nd 

r[i/2] = 0F. 
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